arXiv: 1502.0768lvl [physics.gen-ph] 14 Feb 2015 


International Journal of Modern Physics B. Vol.27. No. 9. (2013) 1330005. 


Review of Some Promising Fractional Physical Models 

Vasily E. Tarasov 

Skobeltsyn Institute of Nuclear Physics, 

Lomonosov Moscow State University, Moscow 119991, Russia 
E-mail: tarasov@theory.sinp.msu.ru 

Abstract 

Fractional dynamics is a field of study in physics and mechanics investigating the 
behavior of objects and systems that are characterized by power-law non-locality, power- 
law long-term memory or fractal properties by using integrations and differentiation of 
non-integer orders, i.e., by methods of the fractional calculus. This paper is a review of 
physical models that look very promising for future development of fractional dynamics. 
We suggest a short introduction to fractional calculus as a theory of integration and differ¬ 
entiation of non-integer order. Some applications of integro-differentiations of fractional 
orders in physics are discussed. Models of discrete systems with memory, lattice with 
long-range inter-particle interaction, dynamics of fractal media are presented. Quantum 
analogs of fractional derivatives and model of open nano-system systems with memory 
are also discussed. 
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1 Introduction 


In this paper we review some applications of fractional calculus and fractional differential equa¬ 
tions in physics and mechanics. The interest in such applications has been growing continually 
during the last years. Fractional calculus is a theory of integrals and derivatives of any arbi¬ 
trary real (or complex) order. It has a long history [1| [2, S] from 30 September 1695, when the 
derivatives of order a = 1/2 has been described by Leibniz in a letter to L’Hospital. Therefore 
this date can be regarded as the birthday of Fractional Calculus. We can probably think that 
Joseph Liouvillc was a first in application of fractional calculus in physics [3]. The fractional 
differentiation and fractional integration go back to many great mathematicians such as Leibniz, 
Liouvillc, Riemann, Abel, Riesz, Weyl. 

All of us are familiar with derivatives and integrals, like first order 

f'{ x ) = D lf ( x ) = ( Jl /)0r) = dx 1 f(x 1 ). 

and the n-th order 

d n 

/<">(*) = D"J(x) = 

nx nxi pXn-l 

(I n f)(x)= / dx i / dx 2 ... / dx n f(x n ), 

Jo Jo Jo 

where n is integer number n = 1,2,..., i.e. n G N. 

Mathematicians consider the non-integer order of the integrals and derivatives from 1695, 

= D°f(x) = ?, (r/)M =?, 


where a G M or a G C. 

At this moments there are international journals such as ’’Fractional Calculus and Applied 
Analysis”, ’’Fractional Differential Calculus”, ’’Communications in Fractional Calculus”, which 
are dedicated entirely to the fractional calculus. 

The first book dedicated specifically to fractional calculus is the book by Oldham and 
Spanier |3j published in 1974. There are two remarkably comprehensive encyclopedic-type 
monographs. The first such monographs is written by Sarnko, Kilbus and Marichev [5j and 
it was published in Russian in 1987 and in English in 1993. In 2006 Kilbas, Srivastava and 
Trujillo published a very important and other remarkable book [6], where one can find a modern 
encyclopedic, detailed and rigorous theory of fractional differential equations. It should be noted 
the books on fractional differential equation by Podlubny |7] and an introduction to fractional 
calculus for physicists by Herrmann [8]. There exist mathematical monographs devoted to 
special questions of fractional calculus. For example, these include the book by McBride [9] 
published in 1979 (see also [IQ]), the work by Kiryakova m of 1993. The fractional integrals 
and potentials are described in the monograph by Rubin [12] , the univalent functions, fractional 
calculus and their applications are described in the volume edited by Srivastava and Owa CBJ. 
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Fractional differentiation inequalities are described in the book by Anastassiou [T4] published 
in 2009. 

The physical applications of fractional calculus to describe complex media and processes 
are considered in the very interesting volume edited by Carpintery and Mainardi [15] publish 
in 1997. Different physical systems are described in the papers of volumes edited by Hilfer 
[15] in 2000, and the edited volume of Sabatier, Agrawal and Tenreiro Machado mi published 
in 2007. The most recent volumes on the subject of application of fractional calculus is the 
volumes edited by Luo and Afraimovich [18] in 2010, and by Klafter, Lim, and Metzler [19] in 
2011. The book by West, Bologna, and Grigolini [20] published in 2003 is devoted to physical 
application of fractional calculus to fractal processes. The first book devoted to application of 
fractional calculus to chaos is the book by Zaslavsky [21] published in 2005. The interesting 
book by Mainardi [22] devoted to application of fractional calculus in dynamics of viscoelastic 
materials. The books dedicated specifically to application of derivatives and integrals with 
non-integer orders in theoretical physics are the remarkable books by Uchaikin [23] 124], and 
the monographs by Tarasov [25] [26]. We also note a new book by Uchaikin and Sibatov [ 27] 
devoted to fractional kinetics in solids. 

Due to the fact that there are many books and reviews on application of fractional calculus to 
describe physical processes and systems, it is almost impossible in this review to cover all areas 
of current research in the held of fractional dynamics. Therefore we must to choose some of the 
areas in this held. We have chosen areas of fractional dynamics that can be considered as the 
most perspective directions of research in my opinion. These areas are not related to a simple 
extension of equations with derivatives of integer order to non-integer. We consider fractional 
models that give relationships between different types of equations describing apparently the 
system and processes of various types. In addition, we think that these areas and models 
can give new prospects for a huge number of fundamentally new results in the construction 
of mathematical methods for the solution of physical problems, and in the description of new 
types of physical processes and systems. 

As a hrst type of models, we consider discrete maps with memory that are equivalent to 
the fractional differential equations of kicked motions. These models are promising since an 
approximation for fractional derivatives of these equations of motion is not used. This fact 
allows us to study the fractional dynamics by computer simulations without approximations. 
It allows us to find and investigate a new type of chaotic motion and a new type of attractors. 

As a second type of promising fractional models, we consider the discrete systems (or me¬ 
dia) with long-range interaction of particles, and continuous limits of these systems such that 
equations of motion with long-range interaction are mapped into continuous medium equations 
with the fractional derivatives. As a result we have microscopic model for fractional dynamics 
of complex media. 

The third type of models is related to the description of the fractional dynamics by mi¬ 
croscopic models of open quantum systems which interacts with its environment. We give an 
example that demonstrate that time fractional dynamics and a fractional differential equation 
of motion can be connected with the interaction between the system and its environment with 
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power-law spectral density. 

We also consider fractional models that allow us to describe specific properties of fractal 
media dynamics; quantum analogs of fractional derivative with respect to coordinate and mo¬ 
mentum; the importance of self-consistent formulation of fractal vector calculus and exterior 
calculus of differential forms that are not yet fully implemented. This review starts with a short 
introduction to the fractional calculus. 


2 Derivatives and integrals of non-integer orders 

There are many different definitions of fractional integrals and derivatives of non-integer orders. 
The most popular definitions are based on the following. 

1) A generalization of Cauchy’s differentiation formula; 

2) A generalization of finite difference; 

3) An application of the Fourier transform. 

We should note that many usual properties of the ordinary derivative D n are not realized 
for fractional derivative operators D a . For example, a product rule, chain rule, semi-group 
property have strongly complicated analogs for the operators D a . 


2.1 A generalization of Cauchy’s differentiation formula 


Let G be an open subset of the complex plane C, and / : G —» C is a holomorphic function. 
Then we have the Cauchy’s differentiation formula 


n\ 


^<*> = 2 ^ 


m 


{z - x) n+l 


dz. 


( 1 ) 


A generalization of ([T]) has been suggested by Sonin (1872) and Letnikov (1872) in the form 


D a J(x) = 


T(a + 1) 

27 ri 


m 


X 


|Q+1 


dz, 


( 2 ) 


where a e R and a ^ —1,-2,— 3,... See Theorem 22.1 in the book by Sarnko, Kilbas, and 
Marichev J5j. Expression (J2]) is also called Nishimoto derivative [28] . More correctly it should 
be called Sonin-Letnikov derivative. 


2.2 A generalization of finite difference 

It is well-known that derivatives of integer orders n can be defined by the finite differences. 
The differentiation of integer order n can be defined by 

£>;/(*) = 


5 






where A 1 / is a finite difference of integer order n that is defined by 




k =0 


( 3 ) 


The difference of non-integer order a > 0 is defined by the infinite series 


- kh). 


(4) 


where the binomial coefficients are 


a 


(T(a + 1) 


p.J r(p + i)r(a-p + iy 

The left- and right-sided Griinwald-Letnikov (1867,1868) derivatives of order a > 0 are 
defined by 


GL r^a 


D a x± f(x ) = Inn 




(5) 


/i-> o h a 

It is interesting that series (J3D can be used for a < 0 and equation (03) defines Griinwald- 
Letnikov fractional integral if 


\f(x)\ < c( 1 + |x|) fi > \a\. 


Then ([5]) can be represented by 


GL r^a 


D a x± f{x) = 


a 


f(x) ~ f{xTz) 


T(1 - a) J 0 

if f(x) G Lp(M), where 1 < p < 1/a and 0 < a < 1. 


yOL-\~ 1 


dz 


2.3 A generalization by Fourier transform 

If we define the Fourier transform operator T by 

1 r+°° 

W )( w)=2 -J f (6) 

then the Fourier transform of derivative of integer order n is 

{FD2f)(u) = {iunrf)(u). 

Therefore we can define the derivative of integer order n by 

D'Jf(x) = 


6 






For f(t) G Li(M), the left- and right-sided Liouville fractional derivatives and integrals can 


be defined (see Theorem 7.1 in [5j and Theorem 2.15 in |6j) by the relations 

{D%f){x) = J- 1 ((±i W )“(JF/)( W )), (7) 

M), (8) 

where 0 < a < 1 and 

( ±iu) a = \u\ a exp(^±sgn(uj)^^J . 

The Liouville fractional integrals (EJ) can be represented by 

(J±/)(x) = —- / z a ~ 1 f(xTz)dz. (9) 

\ a ) Jo 

The Liouville fractional derivatives ([7]) are 

1 d n 

(D±f)(x) = D n x {Il~ a f)(x) = — j^ z^fix =F z)dz, (10) 

where n — [a] + 1. 

We can define the derivative of fractional order a by 

1 r 00 

c D1f(t) = IZ-‘(D!f)(t)=- - - z’'-‘- 1 D"f(x T z)dz, (11) 

i (n — aj j o 


where n — [a] + 1. It is the Caputo derivative of order a [6]. For x G [a, b] the left-sided Caputo 
fractional derivative of order a > 0 is defined by 


c 

a 


D?f(t) 


TTl — CX. 

CL 1 1 


D?m 


1 f l drD^f(r) 
r (n-a)J a (t — r) a ~ n+1 ’ 


( 12 ) 


where n — 1 < a < n, and a /“ is the left-sided Riemann-Liouville fractional integral of order 
a > 0 that is defined by 


= 


f{r)dr 


r(a) J a (■t - r) 


1—a ' 


(a < t). 


Note that the Riemann-Liouville fractional derivative has some notable disadvantages in 
applications such as nonzero of the fractional derivative of constants, 


o 


DfC 


f-ot 

1 r 

r(l-a) ’ 
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which means that dissipation does not vanish for a system in equilibrium. The Caputo fractional 
differentiation of a constant results in zero 


c 

o 


D?C 


0. 


The desire to use the usual initial value problems 


/(h)) = C 0 , o) = C 1 , ( Dlf)(to) = C 2 ,... 


lead to the application Caputo fractional derivatives instead of the Riemann-Liouville deriva¬ 
tive. 

The Riemann-Liouville and Caputo derivative are connected [6|. Let f(t ) be a function for 
which the Caputo derivatives of order a exist together with the Riemann-Liouville derivatives. 
Then these fractional derivatives are connected by the relation 


? = 


*D?m 


m— 1 

E 

k =0 


(t 


\k—a 


T(k — a + 1) 


/“’(«)• 


(13) 


The second term of the right hand side of equation (TT3T) regularizes the Caputo fractional 
derivative to avoid the potentially divergence from singular integration at t — 0. 


2.4 Some unusual properties of fractional derivatives 

Let us demonstrate the unusual properties of derivatives of non-integer orders by using the 
Riemann-Liouville derivatives. 

1) Semi-group property does not hold 


(DZ+DLf)(x) 


1 / 31+1 


wL/x*) - E ( D h‘/)(a) 

k= 1 


(x — a) a k 
T(1 — a — k) 


(14) 


for f(x ) G Li(a,b) and (// + “/)(x) G AC n [a,b], (see equation 2.1.42 in H). 
As a consequence, in general we have 


n a n a v r> 

ly a+ ly a+ W ly a+- 

2) The derivative of the non-zero constant is not equal to zero 

(x — a)~ a 


m + m) = 


T(1 — a) 


(15) 


3) The initial conditions for differential equation with Riemann-Liouville derivative differ 
from the conditions for ordinary differential equations of the integer order 


(oA“ Xr)(0+) = c t , k = 1, ...,n. 


(16) 





For example, conditions (ITfil) for 1 < a < 2 give 

(o Df 1 x)( 0 +) = Ci, (o 2 x)( 0 +) = (o It °A)( 0 +) = c 2 . 

4) Representation in the form of an infinite series of derivatives of integer orders 

r (ct -pi) (x — a) k ~ a 


(d:j)(x) = v 


—' r(fc + l)T(a — k + 1 ) r(fc — a + 1 ) 


Dim 


( 17 ) 


n=0 


for analytic (expandable in a power series on the interval) functions on (a, b ), (see Lemma 15.3 
in [5]). 

5) A generalization of the classical Leibniz rule 


£>“(/») = E 


k =0 


f(n—k)(k) 


from integer n to fractional a contains an inhnite series 

r(a + 1) 


W + (/ 9 ))W = E 


— J T(k + l)T(a — k + 1) 


(Di~*l)(x) Dig 


(18) 


k =o 


for analytic functions on (a, b) (see Theorem 15.1 in [5]). The sum is inhnite and contains 
integrals of fractional order (for k > [ck] + 1). 

6) The increasing complexity of the Newton-Leibniz equation for 

77 / \ Q- fc 

u: + Di + f)(x) = f(x)-Y J r ^y°. +1) 


a+ 


7)(a) 


for f(x) G Li(a,b), (/” a f)(x ) G AC n [a,b ], (see equation 2.1.39 in [6]). 
For 0 < a < 1, we have 

(j“ + z?“ + /)(6) = m - h_AA(/y»/)(«). 


For n G N 


r(a) 

(w/)w = m - E -«)‘ 


k=0 


k\ 


(see equation 2.1.41 in [6]). 

7) For the fractional derivatives there is an analogue of the exponent. The Mittag-Leffler 
function 


EM = E 


T( - ka + 

is invariant with respect to the Riemann-Liouville 

D° + E a [X(x - a) a ] = XE a [X(x - a) c 

(see equation 2.1.57 in [6]). 
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3 Introduction to fractional dynamics 

Fractional dynamics is a field in physics and mechanics, studying the behavior of objects and 
systems that are described by using integrations and differentiation of fractional orders, i.e. 
by methods of fractional calculus. Derivatives and integrals of non-integer orders are used to 
describe objects that can be characterized by the following properties. 

1) A power-law non-locality; 

2) A power-law long-term memory; 

3) A fractal-type property. 

3.1 Fractional diffusion-wave equation 

In mathematics and physics there are well-known the following equations: the diffusion equation 

V 2 w(£, x) = Ci D]u(t, x), 


and the wave equation 

V 2 w(£, x) = C 2 D 2 u(t, x). 

We can consider a generalization of these equations such that it includes derivatives of non¬ 
integer order with respect to time. This generalization describes phenomena that can be char¬ 
acterized by diffusion and waves properties. The fractional diffusion-wave equation is the linear 
fractional differential equation obtained from the classical diffusion or wave equations by replac¬ 
ing the first- or second-order time derivatives by a fractional derivative (in the Caputo sense) 
of order a with 0 < a < 2, 

V 2 w(t, x) = C a q D^u(t, x). 

The solutions of these fractional partial differential equations are described in the book [6] (see 
Section 6.1.2). This equation describes diffusion-wave phenomena [32l [33] . which is also called 
the anomalous diffusion such that we have the superdiffusion for 1 < a < 2, and subdiffusion 
for 0 < a < 1. A more detailed description of these effects and phenomena can be found in the 
reviews [M[55j . 

3.2 Viscoelastic material 

If the force is immediately relaxed, then the deformation disappears. This property is called the 
elasticity. The elasticity is a physical property of materials which return to their original shape 
after they are deformed. The other well-known property is called the viscosity. The viscosity 
of a fluid is a measure of its resistance to gradual deformation by shear stress or tensile stress. 

Mechanically, this behavior is represented with a spring of modulus E , which describes the 
instantaneous elastic response. The stress cr(t) is proportional to the zeroth derivative of strain 
e(t) for elastic solids and to the first derivative of strain for viscous fluids. 
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The elastic solids are described by the Hooke’s law: 

a(t) = E e(t), 


where E is the elastic moduli. 

The viscous fluids are described by the law of Newtonian fluids: 

a(t) = 7 1 D\e(t), 


where 77 is the coefficient of viscosity. 

There are materials that demonstrate these two properties (elasticity and viscosity) at the 
same time. These materials are called viscoelastic. To describe the fractional viscoelasticity G. 
W. Scott Blair (1947) uses the relation |29.13UJ 

cr(f) = E a $Dfe(t), 


where E a is a constant. 

If F(x) is an acting force and x is the displacement, then Hooke’s model of elasticity 

F(x) = —kDf X (t), 

and Newton’s model of a viscous fluid 

F(x) = —kD]x(t), 

can be considered as particular cases of the relation 

F(x) = —kD^x(t). 

This force describe the property that is called fractional friction. 

More complicated fractional models for viscoelasticity of materials are considered in the 
books by Rabotnov [31] and Mainardi [22 . 


3.3 Power-law memory and fractional derivatives 

A physical interpretation of equations with derivatives and integrals of non-integer order with 
respect to time is connected with the memory effects. 

Let us consider the evolution of a dynamical system in which some quantity A(t) is related 
to another quantity B(t) through a memory function M(t) by 



M(t — t)B(t)cIt. 


(19) 


This operation is a particular case of composition products suggested by Vito Volterra. I 11 
mathematics, equation (TT9T) means that the value A[t) is related with B(t) by the convolution 
A(t) = M(t) * B(t). 
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Equation (1T9]) is a typical equation obtained for the systems coupled to an environment, 
where environmental degrees of freedom being averaged. Let us note the memory functions for 
the case of the absence of the memory and the case of power-law memory. 

The absence of the memory: For a system without memory, the time dependence of the 
memory function is 

M(t — t) — M(t) 5(t — t), (20) 

where 5(t — r) is the Dirac delta-function. The absence of the memory means that the function 
A(t) is defined by B(t) at the only instant t. In this case, the system loses all its values of 
quantity except for one: A(t) = M{t)B(t). Using ([T9]) and (120]) . we have 

A(t) = [ M(t)S(t - r)B(r)dr = M(t)B(t). (21) 

Jo 

Expression (l 2 Tj) corresponds to the well-known physical process with complete absence of mem¬ 
ory. This process relates all subsequent values to previous values through the single current 
value at each time t. 

Power-law memory: The power-like memory function is defined by 

M(t — t) = M 0 (t — t) £ ~\ (22) 

where M 0 is a real parameter. It indicates the presence of the fractional derivative or inte¬ 
gral. The integral representation is equivalent to a differential equation of the fractional order. 
Substitution of (| 22 |) into (TT9|) gives the temporal fractional integral of order e: 

A(t) = A I’B(t) = j4) /V - T) c - l B(T)dT, 0 < £ < 1 , (23) 

where A = T(e)M 0 . The parameter A can be regarded as the strength of the perturbation 
induced by the environment of the system. The physical interpretation of the fractional in¬ 
tegration is an existence of a memory effect with power-like memory function. The memory 
determines an interval [0, t] during which B{t ) affects A(t). 

Equation (1231) is a special case of relation for A(t ) and B(t), where A(t ) is directly pro¬ 
portional to B(t). In a more general case, the values A(t) and B(t ) can be related by the 
equation 

f(A(t),M(t) * D™B(t)) — 0, (24) 

where / is a smooth function. For dynamical systems relation (1241) defines a memory effect. 
In this case (|24|) gives the relation f(A(t), q D^B(t)) = 0 with Caputo fractional derivative. 
Relation (|24|) is a fractional differential equation. 
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4 Discrete physical systems with memory 

Discrete maps (universal, Chirikov-Taylor, rotator, Amosov, Zaslavsky, Henon) can be obtained 
from the correspondent equations of motion with a periodic sequence of delta-function-type 
pulses (kicks). 

An approximation for derivatives of these equations is not used. This fact is used to study 
the evolution that is described by differential equations with periodic kicks. 

Example: The universal map without memory 

x n+1 = x n + p n+ iT, p n+1 — p n KT G[x n ] (25) 

is obtained from the differential equation of second order with respect to time 

OO 

D 2 x(t) + KG[x(t)\ S(t/T -k) = 0, (26) 

k =1 

where T = 2n/v is the period, and K is an amplitude of the pulses. 

If G[x] = sin(x), then we have the Chirikov-Taylor map. 

For G[x] = —x we have the Amosov system. 


4.1 Universal map with Riemann-Liouville type memory 

The Cauchy-type problem for the differential equations 

OO 

0 D?x(t) = -KG{x(t)]J25{ T ~k), 1<<*$2 ( 27 ) 

k =1 

with the initial conditions 

(oA“ _1 U(0+) = Cl, (oD?~ 2 x){0+) = ( 0 I?~ a x)(0+) = c 2 , ( 28 ) 

where o-D“ _1 is the Riemann-Liouville derivative, is equivalent to the map equations in the 
form 

x n+i = pk + lVa ( n ~ k + 1 )+ Y(a_ f) ( n + ( 29 ) 

Pn +1 = Pn - KT G[x n ], (1 < a ^ 2 ), ( 30 ) 

where p\ = Ci, and V a (z) = z 0-1 — (z — 1) Q_1 , (z ^ 1). The proof of this statement is gived in 

[39, 253 If G[ x \ — sin(a:), then we have the Chirikov-Taylor map with memory. For G[x\ = —x 

we have the Amosov system with memory. 
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4.2 Universal map with Caputo type memory 

The Cauchy-type problem for the differential equations 

D\x{t)=p{t), (31) 

OO , 

°D?- 1 p(t) = ~KG{x(t)]J2l>U- k )- (1<«< 2 ) (32) 

k= 1 

with the initial conditions x(0) = xo, p( 0) = po, where qD^~ 1 is the Caputo derivative, is 
equivalent to the map equations in the form 


KT a 

x n+1 = x 0 + p 0 (n + 1 )T - —— V' (n + 1 - k) a ~ l G[x k ], (33) 

r <°) hr 

jy^rpa—l n 

Pn+i=Po-f, -ttE[( ,! + 1 - k) a - 2 G[x k ]. (34) 

r <“ - J ) hr 

The proof of this statement is gived in [391 [25] This map allows us to describe fractional maps 
with memory for dynamics with usual initial conditions. 

4.3 Kicked damped rotator map with memory 

Equation of motion 

OO 

0 D?x - q 0 D?x = KG[x\ ^ S(t - nT ), (35) 

n =0 

where gel, 1 < ck ^ 2, /3 = a — 1, and is Riemann-Liouville derivative, can be represented 
in the form of the discrete map 


23+1 


—rr^2Pk+iW a (q,T,n + l - k), 


r(a — i) 


k =0 


Pn+ 1 = e qT (p n + KG[x n ]^, 
where the function W a is defined by 


W a (q, T,m + 1) = T a ~ l [ e~ qTy (m + y) a ~ 2 dy. 

Jo 

The proof of this statement is given in [25] (see also E3BD]. 


(36) 

(37) 


(38) 
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4.4 New type of attractors 

The suggested maps with memory are equivalent to The correspondent fractional kicked dif¬ 
ferential equations [36], 37' j39]. An approximation for fractional derivatives of these equations 
is not used. This fact is used to study the evolution that is described by fractional differential 
equations. Computer simulations of the suggested discrete maps with memory prove that the 
nonlinear dynamical systems, which are described by the equations with fractional derivatives, 
exhibit a new type of chaotic motion and a new type of attractors. For example, the slow con¬ 
verging and slow diverging trajectories, ballistic trajectories, and fractal-like sticky attractors, 
in the chaotic sea can be observed [38] for Chirikov-Taylor map with power-law memory (see 
also mmmmm)- 

5 Dynamics of systems with long-range interaction 

Dynamics with long-range interaction has been the subject of continuing interest in different 
areas of science. The long-range interactions have been studied in discrete systems as well as 
in their continuous analogues. 

The dynamics described by the equations with fractional space derivatives can be charac¬ 
terized by the solutions that have power-like tails [6]. Similar features were observed in the 
lattice models with power-like long-range interactions [45.; 06j 47, 08, [49]. As it was shown 
[5U1 I5T] 1521 [53 , 54; [55, [55] 57[ [58], the equations with fractional derivatives can be directly 
connected to chain and lattice models with long-range interactions. 

Equations of motion of one-dimensional lattice system of interacting particles: 



d 2 u n (t ) 
dt 2 


(39) 


where u n {t) are displacements from the equilibrium, F(u n ) is the external on-site force, and 


OO 




(40) 


n =1 


5.1 Long-range interaction of power-law type 

We define a special type of interparticle interaction 


k^o \k\ a 


dg(k) jq(0) 


A a , a > 0, 0 < \A a \ < oo, 


(41) 


where 



J a (kAx) = ^ e lknAx J(ii) = 2 ^ J{n) cos(fcnAx). 


—iknAx 


(42) 


n=—o o n=1 


n^O 
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This interaction is called the interaction of power-law type a. 

As an example of the power-law type interaction, we can consider 


J{n — m ) = 


| n — m\P +1 


( 43 ) 


The other examples of power-law type interaction are considered in [501 i 5Tl [25] . 

Equations of motion (j39l) with the power-law interaction (l43l) give the following equations 
in the continuous limits: 

1) For 0 < (3 < 2 ((3 ^ 1) we get the Riesz fractional derivative D “ of order a = (3: 


d 2 

-—u(x, t ) — G a A a D^u(x, t) = F ( u(x , t )), 0 < a < 2, (a ^ 1). 

dt z 

2) For f3 > 2 (f3 ^ 3,4, 5,...) we get derivative of order a — 2: 
d 2 


(44) 


dt 2 


u(x, t) + G a ((a — l)D 2 u(x, t) = F (u(x, t )), a > 2, (a ^3,4,...), 


(45) 


where G a = g\ Ax| mm l“ ;2 l is the finite parameter. 

3) For /3 = 1 we get derivative of order a = 1 and ex = 2 for [3 = 3, 5, 7...: 

SHM) iGi a “(f. f ) =F {u{x<t))i 


dt 2 


dx 


(46) 


where G\ = ngAx is the finite parameter. 

4) For (3 — 3, 5, 7... (j3 = 2m — 1, where m — 2, 3,4,...), we get equation with derivative of order 
cx = 2: 

d 2 u(x,t) „ =F{u{xt)): (47) 


where 


G 2 = 


Go 

dt 2 " 9a; 2 

(-l) m - 1 (2vr) 2m - 2 


4 (2m — 2)! B 2 m- 2 g(Ax) 

are the finite parameters, and B 2m -2 are Bernoulli numbers. 

5) For (3 = 2k — 2, where k G N, we have the logarithmic poles. 


(48) 


The effects of synchronization, breather-type and solution-type solutions for the systems 
with nonlocal interaction of power-law type 0 < (3 < 2 (f3 ^ 1) were investigated [53, |5H 
[55. 56]. Nonequilibrium phase transitions in the thermodynamic limit for long-range systems 
are considered in [57]. Statistical mechanics and dynamics of solvable models with long-range 
interactions are discussed in |59j . Stationary states of fractional dynamics of systems with 
long-range interactions are discussed in [58]. Fractional dynamics of systems with long-range 
space interaction and temporal memory is also considered in [55 . [58] . 
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Fractional derivatives with respect to coordinates describe power-law nonlocal properties 
of the distributed system. Therefore the fractional statistical mechanics can be considered 
as special case of the nonlocal statistical mechanics [65]. As shown in the articles [51; [50] the 
spatial fractional derivatives are connected with long-range interparticle interactions. We prove 
that nonlocal alpha-interactions between particles of crystal lattice give continuous medium 
equations with fractional derivatives with respect to coordinates. In the monographs by Vlasov 
[65 ], a nonlocal statistical model of crystal lattice is suggested. Therefore we conclude [25] 
that the nonlocal and fractional statistical mechanics are directly connected with statistical 
dynamics of systems with long-range interactions [ 59] . 


5.2 Nonlocal generalization of the Korteweg-de Vries equation 

The Korteweg-de Vries equation is used in a wide range of physics phenomena, especially those 
exhibiting shock waves, travelling waves and solitons [60]. In the quantum mechanics certain 
physical phenomena can be explained by Korteweg-de Vries models. This equation is used in 
fluid dynamics, aerodynamics, and continuum mechanics as a model for shock wave formation, 
solitons, turbulence, boundary layer behavior and mass transport. 

The continuous limits of the equations of lattice oscillations 


du n (t) 

dt 


+oo 


+ 00 


9 i Ji(n, m)[u 2 n - u 2 n ] + g 3 ^ J 3 (n, m)[u n - u m ], 


m =—oo 
m^n 


m =—oo 
m^n 


give the nonlocal generalization of the Korteweg-de Vries equation 

d d d 3 

—u(x,t) - G 1 u(x,t)—u(x,t) + G 3 —u(x,t ) = 0. 


in the form 


d d ai d a3 

ju(x,t) - G m u(x,t)——u(x,t) + G m -——u(x,t) = 0. 


(49) 


where G ai = gi\/S.x\ ai and G a3 = g 3 |Ax|" 3 are finite parameters. Equation (US]) is a fractional 
generalization of Korteweg-de Vries equation [61] . 

The nonlinear power-law type interactions defined by f(u) = u 2 and f(u) = u — gu 2 for the 
discrete systems are used to derive the Burgers and Boussinesq equations and their fractional 
generalizations in the continuous limit. Note that a special case of this equation is suggested 


m 


5.3 Nonlocal generalization of the Burgers and Boussinesq equa¬ 
tions 

Let us consider examples of quadratic-nonlinear long-range interactions [50; |5T[ [25] . 
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We can consider the discrete systems that are described by the equations 


du n (t) 

dt 


+oo +oo 

= 9i +92 M n -'m)[u n 

m=—oo m =—oo 

m^n rriy^n 



( 50 ) 


where ,J\ (n — m) and J 2 (n — m) define interactions of power-law type with a\ and a 2 . If a.\ — 1 
and CK 2 = 2, then we get the well-known Burgers equation that is a nonlinear partial differential 
equation, which is used to describe boundary layer behavior, shock wave formation, and mass 
transport. If a 2 = a, then we get the fractional Burgers equation that is suggested in [62]■ In 
the general case, the continuous limit gives the fractional Burgers equation in the form 


d_ 

dt 


u(x, t ) + G ai u(x, t ) 


Qai Qa 2 

7T) —; - U(x,t) - G a2 - - U(x : t) = 0. 

d\x\ a ' K ' 2 dx a2 v y 


(51) 


We can consider the chain and lattice equations of the form 


d 2 u n (t) 
dt 2 


+OO +OO 

92 J 2 (n,m)[f(u n ) - f(u m )] +g 4 E 1 / 4 ( 77 ., m)[u n — u m \, 

m =—00 m =—00 

m^n m^n 


(52) 


where f{u) = u — gu 2 , and J 2 (n — m ) and J 4 (n — m ) define the interactions of power-law 
types with ct 2 and cx 4 . If a 2 = 2 and ck 4 = 4, then in continuous limit we obtain the well-known 
nonlinear partial differential equation of forth order that is called the Boussinesq equation. This 
equation was subsequently applied to problems in the percolation of water in porous subsurface 
strata, and it used to describe long waves in shallow water and in the analysis of many other 
physical processes. In the general case, the continuous limit gives the fractional Boussinesq 
equation of the form 

Q2 QC*2 QO-2 Qa 4 

—n(x,t) - G a2 —u{x,t) + gG a2 —u 2 (x,t) + G aA —u( Xl t) = 0. (53) 

Fractional generalization of Korteweg-de Vries, Burgers, Boussinesq equations can be used to 
describe properties of media with nonlocal interaction of particles. 


6 Fractional models of fractal media 

Fractals are measurable metric sets with a non-integer Hausdorff dimension [85l [86]. The main 
property of the fractal is non-integer Hausdorff dimension that should be observed on all scales. 
In real physical objects the fractal structure cannot be observed on all scales but only those for 
which R 0 < R < R rn , where i?o is the characteristic scale of the particles, and R m is the scale 
of objects. Real fractal media can be characterized by the asymptotic form for the relation 
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between the mass M(W ) of a region W of fractal medium, and the radius R containing this 
mass: 

M d (W) = M 0 {R/R 0 ) d , R/Ro > 1. 

The number D is the mass dimension. The parameter D , does not depend on the shape of the 
region W, or on whether the packing of sphere of radius Ro is close packing, a random packing 
or a porous packing with a uniform distribution of holes. 

The fractality of medium means than the mass of fractal homogeneous medium in any region 
W C M n increases more slowly than the n-dimensional volume of this region: 

M d (W) ~ (Vn(kb)) D/n . 

As a result, we can define that fractal medium is a system or medium with non-integer physical 
(mass, charge, particle, ...) dimension. 

To describe fractal media by fractional continuous model, we can use two different notions 
such as density of states c n (D, r) and distribution function p(r). 

(1) The function c n (D, r) is a density of states in the n-dimensional Euclidean space M n . 
The density of states describes how closely packed permitted states of particles in the space 
M”. The expression c n (D,r)dV n represents the number of states (permitted places) between 
V n and V n + dV n . 

(2) The function p(r) is a distribution function for the n-dimensional Euclidean space M n . 
The distribution function describes a distribution of physical values (for example, the mass, 
probability, electric charge, number of particles) on a set of possible states in the space M n . 

Note that some elementary models of fractal density of states and fractal distributions by 
open and closed boxes are suggested in [25] . 

6.1 Homogeneity and fractality 

To describe the fractal medium, we use a continuous medium model. In this model the fractality 
and homogeneity properties can be realized in the following forms: 

(1) Homogeneity: The local density of homogeneous fractal medium can be described by 
the constant density p(r) = p 0 = const. This property means that the equations with constant 
density must describe the homogeneous media, i.e., if p(r) = const and V(W\) = V(W 2 ), then 
M d (Wi) = M d (W 2 ). 

(2) Fractality: The mass of the ball region W of fractal homogeneous medium obeys a power 
law relation M ~ R°, where 0 < D < 3, and R is the radius of the ball. If V n (W\) = \ n V n (W 2 ) 
and p(r,f) = const, then the fractality means that M d {Wi) = \ D M D {W 2 ). 

These two conditions cannot be satisfied if the mass of a medium is described by integral 
of integer order. These conditions can be realized by the fractional equation 

M d (W, t) = [ p(r,t)dV D , dV D = c 3 (D,r)dV 3 , (54) 

Jw 

where r is dimensionless vector variable. 
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6.2 Balance equations for fractal media 

The equation of continuity (mass balance) for fractal media 

(55) 

The equation of momentum balance for fractal media 

P Uk = pfk + V?Pki- (56) 

The equation of energy balance for fractal media 

p(|) e = C 3 1 ( J D,r)c 2 (d, r) p H Viu k + (57) 

Here D is a mass dimension of fractal medium, and 

(J^j = ^ +C 3 1 (D,r)c 2 (d,r)uiV h X7^A = ^(D, r)V fc (c 2 (<i, r)H). 

These equations are proved in [ 1105] . 


6.3 Moment of inertia for fractal bodies 


The moments of inertia of fractal-homogenous rigid ball 


I[°) = 


2D 


r (D) 


-m d r\ 


= 1 + 


2(D-3) 


3(D + 2)‘^" ’ if) ~ ' 3(T> + 2) 
The moments of inertia of fractal-homogenous rigid cylinder 


/(") = 


a 


a + 2 


M d R 2 , 


r(«) 

r(?) 


= 1 + 


a — 2 
a + 2 


(58) 


(59) 


The parameter a: is a fractal mass dimension of the cross-section of cylinder (1 < a ^ 2). This 
parameter can be easy calculated from the experimental data by box counting method for the 
cross-section of the cylinder. 

The periods of oscillation for the Maxwell pendulum with fractal rigid cylinder 


PP = 1 4(q + 1) 

r 0 (2) v 3 (ck + 2) ’ 


(60) 
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The deviation T ( \ (>> from for 1 < a ^ 2 is no more that 6 percent. Equation (I60p allow 
us to use an experimental determination of a fractal dimensional for fractal rigid body by 
measurements of periods of oscillations. 

For a ball with mass M D , radius R. and a mass fractal dimension I), we can consider 
the motion without slipping on an inclined plane with a fixed angle /3 to the horizon. The 
condition of rolling without slipping means that at each time point of the ball regarding the 
plane is stationary and the ball rotates on its axis. The center of mass of a homogeneous 
cylinder moves in a straight line. Using the law of energy conservation, we obtain the equation 

V ( D ) = | 9t sin/3. (61) 

As a result, we have 

f(P) = 21(0 + 2) 

v(3) 5(5£> + 6) [ 1 

Note that the deviation of velocity v(D) of fractal solid sphere from the velocity v(3) of usual 
ball is less than 5 percent. 

The suggested equations allows us to measure experimentally the fractional mass dimensions 
D of fractal materials by measuring the velocities. Note that this measured dimension D must 
be related to the fractal dimension that can be determined by the box counting method. 

6.4 Dipole and quadrupole moment of charged fractal distributions 

The fractional model can be used to describe fractal distribution of charges [87, 88] . The distri¬ 
bution of charged particles is called a homogeneous one if all regions W and W' with the equal 
volumes VoiW) = Vd{W') have the equal total charges on these regions, Qd(W ) = Qd(W'). 
For charged particles that are distributed homogeneously over a fractal with dimension D , 
the electric charge Q satisfies the scaling law Q(R) ~ R D , whereas for a regular n-dimensional 
Euclidean object we have Q(R) ~ R n . This property can be used to measure the fractal dimen¬ 
sion D of fractal distributions of charges. We consider this power-law relation as a definition 
of a fractal charge dimension. If all particles of a distribution are identical, then the charge 
dimension is equal to the mass dimension. In general, these dimensions can be considered as 
different characteristics of fractal distribution. 

Let us consider the example of electric dipole moment for the homogeneous (p(r) = po) 
fractal distribution of electric charges in the parallelepiped region 

0 < x < A, 0 < y < B, 0 < z < C. (63) 

In the case of Riemann-Liouville fractional integral, we have the dipole moment in the 
form ^ b c 

pi D) = [ dx [ dy [ dz x a y a ~ 1 z a ~ 1 = 

r (a) Jo Jo Jo 
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Po(ABC) a A 
T 3 (a) a 2 (a + l)’ 

where a = D/3. The electric charge of parallelepiped region is defined by 


( 64 ) 


Po(ABC)" 

Q(W) - r L D - a?T*(a) ■ 

Then the dipole moment for fractal distribution in parallelepiped is 


pI d) = 




Q(W)A, 


By analogy with this equation, 


P ( „ D) = -4t Q(W)B, pf> = -4 t 0(W)C. 
y a +1 ~ 


ci T 1 


Using a/(a + 1) — D/(D + 3), we obtain 


^ ^ ^ <3) - 


(65) 


( 66 ) 


(67) 


( 68 ) 


where j7 3) = |p^| are the dipole moment for the usual 3-dimensional homogeneous distribution. 
For example, the relation 2 < D < 3 leads us to the inequality 


0.8 < p (D) /P (3) < 1- 


(69) 


These inequalities describe dipole moment of fractal distribution of charged particles in the 
parallelepiped region. 

The example of electric quadrupole moment for the homogeneous (p(r) = p 0 ) fractal distri¬ 
bution in the ellipsoid region W: 


x 2 

A 2 


y 2 z 2 

— H-< 1 

B 2 C 2 ~ 


(70) 


is considered in [25, ETJ [88]. The fractional model of fractal media gives [87j, 88] the electric 
quadrupole moments for fractal ellipsoid 


n (D '> = __n( 3 ) 

^ kk 3D ’ 

n m _ 5tt r 2 (g/6 + i/ 2 ) ^ (3 ) 

v « 0+2 r 2 (£>/ 6 ) v " ’ 

where k y l amd k,l — 1,2,3. For 2 < D < 3, we get 


(71) 

(72) 


l < e27e2 < i, ( 73 ) 

0.6972 < Q<?>/9h < 1. (74) 

These inequalities describe values of the diagonal and non-diagonal elements of the electric 
quadrupole moments for fractal distribution of charged particles in ellipsoid region. 
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6.5 Some applications of fractional models of fractal media 

In this section, we considered some fractional models to describe dynamics of fractal media. In 
general, the fractal medium cannot be considered as a continuous medium. There are points 
and domains that are not filled of particles. We consider the fractal media as special continuous 
media. We use the procedure of replacement of the medium with fractal mass dimension by 
some continuous model that uses the fractional integrals. The main notions that allows us 
to describe fractal media are a density of states and density of distributions. The fractional 
integrals are used to take into account the fractality of the media. Note that fractional integrals 
can be considered as integrals over the space with fractional dimension up to numerical factor 

m 

The suggested fractional models of fractal media can have a wide application. This is due 
in part to the relatively small numbers of parameters that define a fractal medium of great 
complexity and rich structure. In many cases, the real fractal structure of matter can be 
disregarded and we can describe the medium by a fractional model, in which the fractional 
integration is used. The order of fractional integral is equal to the fractal physical dimension 
of the medium. 

The fractional continuous model allows us to describe dynamics of fractal media and fractal 
distributions [ 1481.1251 [26 j. Applications of fractional models to describe fractal distributions 
of charges are considered in [88 ; , :25j. We note that gravitational field of fractal distribution of 
particles and fields can be described by fractional continuous models [106] (see also [81]). Using 
fractional integrals, the fractional generalization of the Chapman-Kolmogorov equation and the 
Fokker-Planck equation for fractal media are derived HE]. We note applications of fractional 
continuous models by Ostoja-Starzewski to the thermoelasticity [89], and the thermomechanics 
[90]. the turbulence of fractal media [91] , the clastic and inelastic media with fractal geometries 
[92], the fractal porous media [93], and the fractal solids [93]. The hydrodynamic accretion in 
fractal media [102 1. [1531 1U3 J is considered by Roy and Ray by using a fractional continuous 
model. 


7 Open system in environment 

The closed, isolated and Hamiltonian systems are idealizations that are not observable and 
therefore do not exist in the real world. As a rule, any system is always embedded in some 
environment and therefore it is never really closed or isolated. Frequently, the relevant envi¬ 
ronment is unobservable or it is unknown in principle. This would render the theory of open, 
non-Hamiltonian and dissipative quantum systems a fundamental generalization of the the¬ 
ory of closed Hamiltonian quantum systems. Now the open, dissipative and non-Hamiltonian 
quantum systems are of strong theoretical interest [ I3Tj - 1561 . 
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7.1 System and environment 

Let Q and P be the self-adjoint operators of coordinate and momentum of the system respec¬ 
tively, and qu and pk describe those of the environment. 

The Hamiltonian H of the system is 

H - = JM + < 75) 

As a model of environment, we consider an infinite set of harmonic oscillators coupled to 
the system. The environment Hamiltonian is 



m nulql \ 

2 )' 


(76) 


This model is called the independent-oscillator model, since the oscillators do not interact with 
each other. 

The interaction between the system and the environment will be considered in the form 


N 

Hi=-Qj2 C nQn 

n= 1 


Q 


c 2 


N 

2 V_— 


(77) 


where C n are the coupling constants. 

Note that the total Hamiltonian H = H s + H e + Hi for the case V(Q) = (MVL 2 /2)Q 2 is the 
well-known Caldeira-Leggett Hamiltonian [67l I681J. 


7.2 Equations of motion for open systems 

Using total Hamiltonian H = H s + H e + Hi , we can derive Heisenberg equations for the system 
and the environment. For the system we have 


dQ 


dP 


N 


1 = aP = M ~ lp • = ~ v 'W + E (' c ' 


dt ih 


n= 1 




y nHn o' 


Q • (78) 


The Heisenberg equations for the environment are 


H = L| q n ,H] =m„Vn, = Pfn.H] = -m n u 2 n q n +C n Q. 


(79) 


Eliminating the operators P and p n , n — 1,..., N, we can write equations (1751) and (170)1 in the 
form 


M 


d 2 Q 
dt 2 


N 


v(Q )=y (c, 


n =1 


a 2 


%qn 


m n u 2 


Q , 


(80) 


24 










(81) 


d 2 1n , 2 _ n n 


The solution of operator equation ( 1ST]) has the form 
q n {t) = q n { 0) cos (u n t) + sin (a : n t) + 


U7/ nfjJrn 


m n (jj. 


Q(t ) sina; n (£ — r) dr, 


(82) 


n^n J 0 


where g n (0) and p n (0) are the initial values of coordinate and momentum operators of the 
environment n-th oscillators. 

Using solution of (l82|) we can derive the equation 


(83) 


+ [ M(t- t + V(Q) = F(t), 


dt 2 


dr 


where the function 


N 


M{t) = Y J 


Cl 


71=1 


■ cos(o; n f) 


(84) 


is called the memory kernel. The one-parameter operator function 


F(t) = V (c n g n (0) cos(ui n t) + sin(o; n f) 


n=l 


n 2 

n Q(0) cos(u n t) 


m n ul 


(85) 


can be interpreted as a stochastic force since the initial states of the environment are uncertain 
and it can be determined by a distribution of the average values of g n (0) and p n ( 0). 


7.3 Quantum dynamics with memory 

The memory function Ai(t) describes dissipation if A 4(t) is positive dehnite and decreases 
monotonically. These conditions are achieved if N —>■ oo and if C1 /and uj n are suffi¬ 
ciently smooth functions of the index n. 

For N —* oo, the sum in equation (|84|) is replaced by the integral 

A4(t) — — [ ^ ^ cos(ujt)duj , ( 86 ) 

^ J—oo ^ 

where J(oj) is spectral density. We assume that the oscillator environment contains an infinite 
number of oscillators with a continuous spectrum. 

For the spectral density 


J(u>) 


N 

7 TU \ -> 


Cl 




5(uj — uj n ), 


(87) 
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equation (| 86 l) gives the memory function (j84l) . If we consider the Cauchy distribution J(u>) = 
a/(co 2 + A 2 ), then equation (1861) gives the exponential memory kernel M.{t) = (a/A) e~ xt . 

We can consider a power-law for the spectral density: 

J{u) = AsJ, 0 < 0 < 1, (88) 


where A > 0. Note that density (I 88 |) leads to the power-law for the memory function Ai(t) ~ 
i _/3 . Equation (188|) can be achieved by a different type of combinations of coupling coefficients 
C(lo) and density of states g(oj ) 

JM = ™g(u)C(u). 

Using the Fourier cosine-transform 


we get the equation 


x a cos (xy)dx = 
r d 2 Q , 


7T 


M- 


2T(a) cos(7tq:/2) 

A 


if , (o < a < 1), 


o D?Q + V'(Q) = F(t), 


dt 2 sin(7T/3/2) 
where 0 -Df is the Caputo fractional derivative 

1 f (. D l Q){r)dr 


o D?Q(t) = 


r(i - p) Jo (t - t )p 


, (o < /3 < 1). 


(89) 


(90) 


(91) 


As a simple example of quantum system, which is described by (IHUlh we can consider the 
linear fractional oscillator that is an object of numerous investigations [108 . 109 . 110111111,11121 
111311114111151111611117111181111911120] because of different applications. 

If F(t) 7 ^ 0, then equation (190T) can be considered as a fractional Langevin equation. We 
also can consider a general quantum analogs of fractional Langevin equation [ 121 . 122; 1123[ 
1124111251112611127111281112911130j that can be connected with the quantum Brownian motion 
that is considered by Lindblad [138]. 

As a result, we obtain a fractional differential equation for operators Q(t) from the inter¬ 
action between the system and the environment with power-law spectral density [ 66 ]. The 
parameter a can be used to control quantum dynamics of nano-systems like individual atoms 
and molecules in an environment [zatsuzacE]- Quantum control is concerned with active 
manipulation of physical and chemical processes on the atomic and molecular scale. Controlled 
manipulation by atomic and molecular quantum systems has attracted a lot of research at¬ 
tention in recent years [69l[70l[n]. Note that the models to control of open quantum system 
dynamics is a very important subject of nanotechnology. 
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7.4 Quantum analogs of fractional derivatives 

One of the way to derive a quantum description of physicsl systems is an application of a 
procedure of quantization to classical models. We can use the Weyl quantization to obtain 
quantum analogs of differential operator of non-integer orders with respect to coordinates. 

The Weyl quantization tt w is defined by 

7T W (<lkA(q, p)) = ^ (QkA + AQf^j , tt w ( p k A(q , p) ) = i (p k A + AP k ), (92) 

*w(Dl k A(q,p)) = -^(PkA-APk'j, ii w [p l pk A{q,p)^ = ^(Q k A- AQ k ^j, ( 93 ) 

for any A = A(Q,P) = 7r w (A(q,p)), where k = 1, Q k = ir w (qk), and P k = n w (p k ). 

Weyl quantization ir w maps [S¥j the differential operator C[q,p, D^, D l p \ on the function 
space and the superoperator C[Lq, Lp, —Lp, Lq) acting on the operator space, where 

L~ a B = —(AB — BA), L+B = -{AB + BA). 
xTi 2 

If A(x) is an analytic function for x G (0, b ), then the Riemann-Liouville fractional derivative 
can be represented in the form 


0 D^A(x) = ^ a(n, a)x n a D™A(x), 

n =0 


(94) 


where 

, \ _ _ r(o? + 1) _ 

r(n + l)T(a — n + l)T(n — a + 1) 

Equation (jUTj) defines a fractional derivative on operator algebra [ BTj . 

The Weyl quantization of the Riemann-Liouville fractional derivatives with respect to phase- 
space coordinates [H] gives 

OO 

„P^ = 7r( 0 D«) (95) 

71=0 

OO 

= -k ( 0 D« ) = £>(«,<*) (L+J""« (L 5 J”, (96) 

71=0 

where L A are dehned by the equations 

L a B = -^{AB-BA), l+b = ^(Ab + bA). 
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For example, we have 


o 


V a Q Q n 


+ X ) n n-a 

r(n + l-ay ’ 


0 


v ap r 


r (^+!) pn- a 

r(n + 1 — a) 


where n ^ 1, and a ^ 0. 

It allows not only consistently formulate the evolution of such quantum systems, but also 
to consider the dynamics of a wide class of quantum systems, such as the non-local non- 
Hamiltonian, dissipative, and nonlinear systems. Quantum analogs of the non-local systems 
with regular and strange attractors can be described |84j . The regular quantum attractors 
can be considered as stationary states of non-Hamiltonian quantum systems. The condition 
given by Davies [ 137 ] defines the stationary state of non-Hamiltonian quantum system. An 
example, where the stationary state is unique and approached by all states for long times is 
considered by Lindblad j!38j for Brownian motion of quantum harmonic oscillator. Spohn 
[ 1391, HO, ljl] derives sufficient condition for the existence of a unique stationary state for the 
non-Hamiltonian quantum system described by Lindblad equation. The stationary solution of 
the Wigner function evolution equation for non-Hamiltonian quantum system was discussed in 
[H 2j, US]. Quantum effects in the steady states of the dissipative map are considered in |144j . 
Stationary pure states of quantum non-Hamiltonian systems are considered in [ 1451 146 j. For 
classical non-Hamiltonian systems, stationary states are presented in [ 1481447 :, ll49j . 


8 Fractional generalization of vector calculus 

The fractional calculus has a long history from 1695. The history of fractional vector calculus 
is not so long. It has little more than ten years and can be reduced to small number of papers 
(about twenty articles) (see Refs, in [64]). 

A consistent fractional vector calculus is important for application in the following research 
directions. 

1) Non-local statistical mechanics. It should be noted book by Vlasov [65] book is entirely 
devoted to nonlocal statistical mechanics. The fractional derivatives in equations can be con¬ 
nected with a long-range power-law interparticle interactions in statistical mechanics [591150]. 

2) Non-local electrodynamics [76] [77] and [78], where the spatial dispersion describes non¬ 
local properties of media. 

3) Non-local hydrodynamics and waves propagation in media with long-range interaction 

[79] (see also Sec 8.16 in [25]). 

It is known that the theory of differential forms is very important in mathematics and physics 

[80] . The fractional differential forms can be interesting to formulate fractional generalizations 
of differential geometry, including symplectic, Kahler, Riemann and affine-metric geometries. 
These generalizations allow us to derive new rigorous results in modern theoretical physics and 
astrophysics [81], and in fractional generalization of relativistic field theory [ 150 .. 15111152111531 
11541 1155111561 11571 115811159111601 11611 116211163] in curved space-time. We assume that the 
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fractional differential forms and fractional integral theorems for these forms can also be used 
to describe classical dynamics [82] and thermodynamics. 

It is important to have fractional generalizations of the symplectic geometry, Lie and Poisson 
algebras, the concept of derivation on operator algebras. It allows to apply a generalization of 
algebraic structures of fractional calculus to classical and quantum mechanics. Note that the 
theory of operator algebras is very important in quantum theory 183] . 
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